HIGHER REGULARITY OF HOLDER CONTINUOUS SOLUTIONS OF PARABOLIC EQUATIONS 

WITH SINGULAR DRIFT VELOCITIES 



SUSAN FRIEDLANDER AND VLAD VICOL 

ABSTRACT. Motivated by an equation arising in magnetohydrodynamics, we prove that Holder continuous weak solutions 
of a nonlinear parabolic equation with singular drift velocity are classical solutions. The result is proved using the space-time 
Besov spaces introduced by Chemin and Lerner, combined with energy estimates, without any minimality assumption on the 
Holder exponent of the weak solutions. 



1. Introduction 

In this paper we address the smoothness of Holder continuous weak solutions of the scalar nonlinear parabolic 
equation with singular drift velocity 

d t 6 - A6 + (u ■ V)0 = (1.1) 

uj = (1.2) 

on R d x (0, oo), where d>2, {Ty^f . =1 is a d x d matrix of Calderon-Zygmund singular integral operators, and the 
summation convention on repeated indices is used throughout. The drift velocity is taken to be divergence-free, i.e. 

V-u = (1.3) 

which is ensured by (11 .2b if the matrix {T^} is taken such that didjTijf = for smooth functions /. The system 
( 11.1 b — dTT3T > is supplemented with the initial condition 

0(-,O) = o (1.4) 

where O £ L 2 (K d )has zero mean on M. d . We note that J Rd 9(x, t) dx is conserved in time by solutions of dl.U - fTzl . 

The global in time existence of finite energy weak solutions to (11.1 b — ( TT~4-b has been proven by the authors of this 
paper in Ifl2l . Additionally, in IPT21 we prove that for positive time the weak solutions are in fact Holder continuous 
(see also [|20| ). In the present paper we address the higher regularity of these Holder continuous weak solutions, 
by proving that they are classical solutions (even C°° smooth) for positive time. This result was announced in ifTZl 
Lemma 3.4]. 

The motivation for studying advection-diffusion equations with drift velocities as singular as those given by ( 11.21 
came from the three-dimensional equations of magneto-geostrophic dynamics, for a rapidly rotating, electrically con- 
ducting fluid (cf. Moffatt |[T8l ). A well studied advection-diffusion equation also arising in rotating fluids (cf. Con- 
stantin, Majda, and Tabak (9)) is the so-called critical surface quasi-geostrophic (SQG) equation. This 2-dimensional 
equation has the form 

a t 6» + (-A) 1/2 6» + (ii-V)c? = (1.5) 
u = V x (-A)- 1 ' 2 0. (1.6) 

Although there are significant differences between the systems fll.l| i-( 11.4t and ( |1.5t -( [L6] l, in both cases L°°(M d ) is 
the critical Lebesgue space with respect to the natural scaling of the equations. The criticality of the L°° norm with 
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respect to scaling also holds for the modified critical surface quasi-geostrophic equation considered by Constantin, 
Iyer, and Wu in JTU) 

d t 9 + (-A)^ 2 9 + (u • V)0 = (1.7) 
u = \7 ± (-A} 3 / 2 - 1 e (1.8) 

where (3 G (0, 1). In a recent paper, Caffarelli and Vasseur |JT] used De Giorgi iteration to prove that weak solutions 
of dl.5t -( [TT6l i, with L 2 initial data, are smooth for positive time. A different proof of global regularity for (11.5t - dT3t 
was given independently by Kiselev, Nazarov, and Volberg lfl3l (see also Kiselev and Nazarov fT4l ). The proof of 
Holder regularity of weak solutions to (ll.lt - dl.4l given by the authors of the present paper in llTZl is also based on 
a suitable modification of the De Giorgi method, along the lines of [1|. Once the weak solutions to (ll.lt - dl.4t are 
Holder continuous, we expect to be able to bootstrap to higher regularity, since the Holder C a norm is subcritical with 
respect to the natural scaling of the equations, for any a > 0. This matter is however not automatic due to the singular 
velocity drift u, which by (TTT2T > lies in LfC^ 1 H L 2 x , whenever 9 G LfC^ n if i? The following theorem is the 
main result of the present paper. 

Theorem 1. Let 9 Q G L 2 (R d ) be given. Let 

9 G L°°([0, oo); L 2 (R d )) n L 2 ((0, oo); H 1 ^)) n L°°((0, oo); C a (R d )) (1.9) 

be a Holder continuous weak solution of d 1 . 1 t -( TT~4t , evolving from 9q, where a G (0, 1) is given. Then, the solution 
is classical, that is 

0eL oo ([t o> oo);C 1 ' J (R ,i )) (1.10) 

for any to > and 6 G (0, 1). 

The issue of proving higher regularity of Holder continuous solutions to a fractional advection-diffusion equation 
has been considered in the context of the supercritical SQG equation 1111 . the critical SQG equation (fl~), the modified 
critical SQG equation ifTOlfTTl , and in the recent preprint lED which addresses a linear equation with singular drift. The 
natural characterization of Holder spaces in terms of Besov spaces were utilized in ifTOl for dl.7t - dl.8t . respectively in 
ifm for the supercritical SQG equation, to prove that if a solution is C a for some a G (0, 1), then in fact it lies in a 
more regular Holder space, which can be bootstrapped to prove the classical solution is classical. 

The techniques used in ifTUl [TT1 may be applied in order to prove higher regularity for the system (ll.lt - (11.4t . but 
only once the C a regularity is such that a > 1/2 (this was also pointed out in ifTTl ). However, if we only know that 
a weak solution of dl.lt - dl.4t is in C a with a G (0, 1/2), the velocity field is too rough, and the aforementioned 
method of ITOl ITT1 does not apply directly. We find that it is necessary to use different arguments to obtain the desired 
result. By working in the Chemin-Lerner space-time Besov spaces (see JSJ) we make use of the smoothing effect of the 
Laplacian at the level of each frequency shell, which enables us to take advantage of the extra a priori information that 
u G L\ x . The principal difficulty lies in treating the high-high frequency interaction in the paraproduct decomposition 
of the nonlinear term (cf. d4.5t below). The main result of this paper is the proof of higher regularity for solutions 
of dl.lt - dl.4t without any minimality requirement on a > 0, and in any dimension d > 2 (cf. Theorem [U. The 
method introduced in order to prove Theorem [T]also gives new higher regularity results for the system (11.7t - (11.8t , in 
the parameter range j3 G (1, 2) (cf. Theorem[6]below). 

This paper is organized as follows. In Section|2]we recall some facts about Besov spaces. Section[3]gives the proof 
of Theorem Q] for a > 1/2, while for the case < a < 1/2 the proof is given in Section |4] Section [5] contains a 
description of our results for the modified critically dissipative SQG equations. 

2. Preliminaries 

Let {(f>j}j<zz be a standard dyadic decomposition of the frequency space R d , with the Fourier support of the 
Schwartz function fa being {2 3 " 1 < |£| < 2 3 + 1 }, and where = 1 on R d \ {0}. As usual, define 

Ajf = fa * f and Sj = J2k<j ^jf f° r a ^ Schwartz functions /. 
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For s G R and 1 < p, q < oo the homogeneous Besov norm B & p is classically defined as 

ll/ll^ tg HI 2J WlMI w (2.D 

whenever q G [1, oo), while in the case q = oo one defines 

=su P 2^||A j /|| L p. (2.2) 

Recall that L°° n B^ x = C s is the Holder space with index s, except when s is a nonnegative integer (then we 
recover the Zygmund spaces C*). For any r G [1, oo] we classically let L r (T; B p ) denote the set of all Bochner 
i r -integrable functions on I, with values in B p , where I C [0, oo) is some given time interval. 

Lastly, for s G R, a time interval I, and 1 < r, p, q < oo we recall the Chemin-Lerner space-time Besov spaces 
L r (I; B^ ), with norm 

(2.3) 

«(Z) 

with the usual convention of taking a supremum in j if = oo. Note that L r (l; B pr ) = L r (I; B p r ) for all r > 1. 



2 JS 



l|A 3 -/(-,*)lli- 



l/r 



3. Proof of the main theorem in the case a e (1/2, 1) 

In this case the proof follows directly from iTTOl [TTI . with only slight modifications, so we give very few details. 
First, note that if 9 is as in the statement of the lemma, then 9 G L°°([to, oo); Bp^o), where a p = (1 — 2/p)a, and 
p G [2, oo) is fixed, to be chosen later. Then, for j G Z fixed, we have 

— HA^Hk + J \A j 0\*-*A j 0(-A)A j 0= -J \A j 9\p- 2 A j 9A 3 (u-V9). (3.1) 

Upon integration by parts and using |[T6l Proposition 29.1] (cf. |[T9l , see also J7]|23l for the fractionally diffusive case), 
the dissipative term is bounded from below as 

/o 2 i 
|A^r 2 A^(-A)A^ dx > — ||A^||^, (3.2) 

where C = C(d,p) > is a sufficiently large constant. The main difficulty lies in estimating the convection term. 
This is achieved in ifTOlfTTl by using the Bony paraproduct decomposition 

Aj(u-V6)= J2 A i V ■ (Sk-iuA k 6) + J2 A j( A kU ■ VS fc _i0) + J2 J2 A,V-(A fe uA^). (3.3) 

|j-fe|<2 |j-fc|<2 k>j-l |fc-{|<2 

When integrated against Aj6\ Aj9\ p ~ 2 , (|3.31 l gives rise to three terms on the right side of (13.1) . The first two terms 
(when \j — k\ < 2) may be bounded favorably for any a > 0, by first integrating by parts the derivative contained in 
Sk-iik = djSk-iTjj9, then using a commutator estimate, the Holder and Bernstein inequalities (see ifTUl for details). 
However, the third term on the right side of (13.3) gives rise to an integral which may only be bounded favorably when 
a G (1/2,1). Indeed, from the Holder inequality we obtain 

< C\\A 3 9\\ p - p 2 V J2 E IIAHMIA^IUoo 
k>j-l |k-J|<2 

< C7||A^||^ 2 2-''||0|| CQ 2 fe(1 - a - a ^2 te *||A fe 0|| L p. (3.4) 

k>j-i 

Given that 9 G Bp^o, the sum of the right side of the above estimate is finite only if a + a p > 1. The later holds if and 
only if a > 1/2 (and p is large enough, depending on a). However, if a G (0, 1/2) it seems that the method of ifTOl 



E E 

I, ->» A 1 I 7_ 1 I *s n 



(A fe uA^)A^|A^| 



p-2 
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cannot be applied directly. We overcome this difficulty in Section |4]below. In the case a £ (1/2, 1), the right side of 
(13. 4t does remain bounded and the estimate on the nonlinear term may be summarized as 



Combining (13. lb — (13.5b with the Gronwall inequality, and taking the supremum in j, we obtain that 

9 £ L°°([ti, oo); Bp a £(R d )) 



(3.5) 



(3.6) 



foranyii > to- Using the Besov embedding theorem we obtain that 9 £ L°°([ti, oo); B^,oo P (K d )), for any t\ > to, 
where e p = (4a + d)/p < (4 + d)/p. Letting p > (4 + d)/(2a — 1) concludes the proof of the theorem in the case 
a £ (1/2,1). 

4. Proof of the main theorem in the case a £ (0, 1/2] 

Let us fix X = [to, t{\, for some < to <t\. The following lemma gives the principal estimate needed in the proof 
of TheoremQ] 



Lemma 2. Let 9 be a weak solution of (ll.l| i-( fl~4] i which is Holder continuous, that is 

9 £ L 00 (l;L 2 (m d ))nL 2 (l;H 1 (R d ))nL 00 {l; C a (R d j) 
for some a £ (0, 1/2). If additionally 

9eL 2 (l;Bl 2 (R d )), 

for some p > 2, then we have 

9eL 2 (l; Bl r (R d )) 
for all 1 < r < oo, and for all q £ (p, m a p), where m a = (1 — a)/(l — 2a) > 1. 



(4.1) 
(4.2) 
(4.3) 



Proo/ Apply Aj to ( fTH . multiply by Aj0| Aj0| 9-2 , integrate over R d , and use HU Proposition 29. 1] (cf. ATI [1911231 ). 
to obtain 



~||A J (9||i 9 + C 2^-||A J fl||«,< 



Aj(u- V0)A j 0\A j 9\ q - 2 



(4.4) 



for some constant c = c(d, q) > 0. Using the Bony paraproduct decomposition and the divergence free nature of u 
(and hence of Sk-iu and Ajt-u) we write 



A J {u-V9)= A,-V • (5 fc _iuA fc fl) + ]T A i (A fc «-V5 fc _ifl)+ ]T ]T A,-V • (A fc uA, 

b'-fc|<2 |j-fc|<2 k>j-l |fc-{|<2 



From the Holder inequality, ( |4.4t , and (14.5) we hence obtain 

1 d 



qdt 



A 3 -0||£, + c2 2 ^||A^]|| 5 < CIIA^H^ 1 (Ji + J 2 + J 3 ) 



where c = c(<i, p) > is a sufficiently small constant, and we have denoted 

Jx= J2 l|A,-V-(5 fe _iuA fe e)||^ 

U-fc|<2 

J 2 = ^ ||Aj(A fc u-VS fc -ifl)||L« 

|j-fe|<2 

J 3 = ^ ^ ||A 3 V-(A fc wA 



(4.5) 

(4.6) 

(4.7) 
(4.8) 
(4.9) 



fc>i-i |fe-/|<2 
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We bound J x using the Bernstein inequality, the boundedness of Calderon-Zygmund operators on L q , and the triangle 
inequality < Y^i<k-i \\^i u \\li for all q S [1, oo] (note that u(0) = 0), to obtain 

,h<cv J2 E ii^i«iu«ii^iu- 

\j-k\<2 Kk-l 

<CV J2 2- to (2 feQ ||A^||^) 2 '||A^||l, 

|j-fc|<2 Kk-l 

<C\\e\\ c ^ Yl 2 ~ ka E 2'(2 / ||A^|| LP ) p/9 (2 /Q ||A^|| i o C ) 1 " p/9 2-'^ +a ( 1 -^» 

\j-k\<2 Kk-l 

< C\\ef c - p/q 2 j ^ 2- ka 2 lil -P/«- a(1 -P/i» (2 l \\A l 6\\ LP ) p/q 

\j-k\<2 Kk-l 

< c\\e\\ 2 jl /q 2^~ a ^ 2l(1 ~ p,q ~ a(1 ~ p/q)) (2'l|A ; e|kO P/9 - ( 4 - 10 ) 

In the above estimate we also used the interpolation inequality 1 1 / 1 1 < ||/|| p / q \\f\\ ]^'^ q , which holds for all functions 
/ G L p n L°°, and any q 6 (p, oo). Note that since a < 1 we have 



l---a(l--)>Q<*l-->0o-q>p, (4.11) 



so that for all q > p and all s £ [1, oo) we have 



l/s 

2sl{l-p/q-a{l-p/q)) 1 < (J2i( 1 -p/l- a { 1 -p/l)) _ (4.12) 



We bound J 2 similarly, 



J2 ^ C £ E l|Afc«||£«||VA,fl||io. 

|j-fe|<2 i<fe-l 

<C ^ 2 fe ||A fc 0|| L9 £ 2'HA^IUoo 

|j-fe|<2 Kfc-1 

<C ^ 2 fc ||A fe 0||^||A fc 0||i- p/9 Y, 2l[1 ~ a) (^IIA^IU-) 

\j-k\<2 Kk-l 

< C\\6\\ c - E 2 * (2 fc ||A fe 0|| LP ) p/9 (2 ka \\A k e\\ L ^ 1 ~ p/q 2- k ^ + ^ 1 ~P/^2 k ^ 

\j-k\<2 

< c\\6fcJ /q 2 j( - 2 - a - p /< 1 - a ( 1 - p / q » {2 k \\^kO\\ LP ) P/q ■ (4.13) 

|j-fe|<2 
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Note that here we used a < 1 to obtain that J2i <k -i 2 /(1 ~ q) < C2 ki - 1 ~ a \ Lastly, we bound J 3 as 

J 3 <C2 j l|A fc u|| is ||A^||i=c 

k>j-l |fc-Z|<2 

<CV 2 fc ||A^|| L9 2- iQ (2' Q ||A^|| L oo) 

k>j-l |fc-(|<2 

<CV\\0\\c~ £ 2 fe ( 1 - Q )||A^||^||A^||i- p/9 
k>j-i 

< CV\\0\\ c ~ 2M1 " a) (2 fc ||A fc 0|| LP ) p/<? (2 fcQ ||A fc 0|| L o.) 1 " p/<? 2- fc ^+ a ( 1 -P/ < ?)) 

fc>j-i 

< C2m9\\ 2 c - a p/q 2 k{1 ~ a ~ p/q ~ a{1 ~ p/q)) (2 k \\&kO\\L P y /9 ■ (4.14) 

k>j-l 

Here as before we used the Bernstein's inequality and the bound || Afcu||ig < C2 k \\A k 0\\L<i- If we let 

1 - a 

p<q<m a p=- — p, (4.15) 

1 — 2a 

the exponent of 2 k in the last inequality of ( 14.141 ) lies in the range 

-a<l-a---a(l--) < 1 - — - a ( 2 ) = (4.16) 

q \ qj m a \ m a J 

since a G (0, 1/2). Therefore, if q G (p, m a p), for any s G [1, oo) we have 

2 sfc (i-"-p/9-"(i-p/9)) J < (72 J ( 1_Q_p / 9_Q ( 1_p / |? ^ (4 17) 

K k>j-i J 
We insert the bounds ( 14.10b . ( 14.13b , and (14.14b into ( 14.6b and obtain the a priori estimate 

||^6>|| ig + c2 2 ^ H^^^IU. < Cll^H^-^/^a^f 1 --) s^^Ci-^/^-^Ci-^/^)) (2 k \\A k 6\\ LP ) p/q 

+ c\\e\\ 2 c l /q v l - 2 ~ a -p'^ 1 'Pi^ J2 (2 fc ||A fc #|| LP ) p/<? 

|i-fc|<2 

+ C\\0\\ 2 Ca p/q 2 : > 2 k ^- a - p ^- a< - 1 -^^ (2 k \\A k 9\\ LP ) p/q . (4.18) 
k>j-i 

We apply Gronwall's inequality and obtain 

IIA^^IU^e-^^^HA^^o)!!^ 

+ ^ll e ll^ P (iV) 2J( ^ Q) E 2fe(1 ^ / ^ Q(1 " P/9)) f e- c(t - s)223 ( 2 k \\A k 9(s)\\ LP ) p/q ds 

k<j 

+ ^ll^ll^ P (i;C°) 2 ' (2 " Q " P/9 " Q(1 " P/9)) E f e~ c{t - s)22j (2 k \\A k 8(s)\\ LP ) p/q ds 

+ C\\e\\ 2 -»l q . ca) V 2 k ^- a -P^ a ^P/^ f* e -c(t-s)2V (2 k \\A k 9( s )\\ LP ) p/q ds. (4.19) 

k>j-l *<> 

Using the Young-type inequality 

\\f*9\\L* { T) < \\fhHT)\\g\\L>m ^ 11/11™ llsllw^m^ 29 (4.20) 
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and the bound 



with r = 1, we obtain that 

rt 



-c2 2 ^(i-t ) 



Lr{X) 



< min{C2- 2j /MZ| 1 / r } 



(4.21) 



f e _c(t-. )a « (2^||A^( s )|| iP ) p/9 d S < Cmin{2- 2 ^, (2 k \\A k 9\\ LHX , LP) ) p/q |J|(^)/ 2 «. (4.22) 

"'*o L 2 (I) 

We take the L 2 (I) norm of ( 14.19t , use the bound (I4.22l i above, combined with the discrete Holder inequality, the fact 
that 9 G L 2 (l; B^ 2 ) = L 2 (l; B^ 2 ), the estimates (14.12b and (14.17b with s = 2q/(2q - p), and obtain that 

+ c\\ef-Jll iCa) \\e^ {x ^ (^m-^-p/i-<i-v/i)) m m{C2~ 2j , |x|}) 

< C||^o)||^||^(t )|lc« P/9 (2-^^) min{2-j, XI 1 / 2 }) 

+ ^l^ll^i^ll^ll^x;^ 2) \X\ (q ' p)/2q ( 2 J( 2 - Q -f/9-«( 1 -f/?)) min{C2-^', |X|}) (4.23) 

for all q G (p,m a p). Note that #(to) G £ p since we a priori have 9 G L^°L 2 n LfL™. We multiply the above 
estimate on both sides by 2 J and take an £ r (Z)-norm, to obtain that 

\\ e \\L*{X;Bl r ) = ll 23 ll A J' e lli a (X;i«)ILr( Z ) 

< ll^o)||^||e(t )|lc^ /9 j^^ 1 -^ 1 -^)) min{2^, l^| 1/2 }[ r(z) 
+ C||^||^V)ll^lli / / tt .Bi J27| (9_p)/2g 2 3(3 " Q " p/ ^ Q(1 " p/9)) min{C72- 2 Ml|} . (4.24) 

The key observation is that for all q G (p, m Q p) we have — a(l — p/g) < 0, 1 — a(l —p/q) > 0, 1 — a — p/q — Qf(l — 
p/g)) < 0, and 3 — a — p/q — a(l — p/g)) > 0, so that the two l r norms on the right side of the above estimate are 
finite, for any 1 < r < 00. We have hence proven that 

0eL 2 (l;Bl r ) (4.25) 

for any 1 < r < 00, and any q G (p, m a p), concluding the proof of the lemma. □ 

The following lemma shows how one may bootstrap the arguments in Lemma|2]in order to control 9 in L 2 W^'°°. 

Lemma 3. Let 9 be a weak solution of (ll.lt -( 11.4t which is Holder continuous, that is 

9 G L°°(l; L 2 (R d )) n L 2 (l; S rl (R d )) n L°°(X; C* Q (R d )) (4.26) 

for some a G (0, 1/2). Then we have 

W6 e L 2 (Z;L° (R d )). (4.27) 

Proof. We note that H 1 = B\ 2 so that we may apply Lemma[2]with p = 2, and obtain that 9 G L 2 (I; B^ 2 ) for any 
g G (2, 2m a ). Since m Q > 1, we may bootstrap and apply Lemma[2]once more to obtain that 9 G L 2 (I; B^ 2 ) for all 
g G (2, 2to 2 ). For any fixed p > 2, since m k a diverges as k — » 00, we may iterate Lemma |2] finitely many times and 
obtain that 9 G L 2 (X; r ), for all r G [1, 00]. 

Fix p large enough, to be explicitly chosen later, and let q = p(l + m a )/2. From the estimate ( 14.23b , for any e > 
we have 

2^)||A J -0|| i2(Z;Lg) 

<C\\9(t )\\% q \\9(t Q )\\ 1 -J ,/q min{C2-''(^ Q ( 1 -^),|X| 1 / 2 2^ 1 + e - Q ( 1 -f/9))} 

+ C\\9\\ 2 L J!ll A\9\\ p,q , \l\(v-p)nq m i n {C2^ f - +1 -P/v-^ 2 ~P/^\2^^ (4.28) 
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where q = p(l + m a ) /2. We now pick a suitable e > 0, so that the i 1 (1) norm of the right side of (|4.28t is finite. For 
this to hold, we need that the following four bounds to hold true 

e - a ( 1 - -^—] <0^e< (4.29) 

2 - 3a - a 2 

1 + e - a 1 - > e > — (4.30) 

2 — 6a 



1 


f TO a 




2 


1 


f m Q 




2 


1 


f m Q 




2 



2 (~ 2 \ „ (l-2a)(2-3a-a 2 ) 

m Q \ 1 + m a J (1 — a)(2 — 3a) 

2 (~ 2 \ n 2 - 3a + a 2 -2a 3 
e + 3 a 2 - > ^ e > — — — — (4.32) 

m a \ 1 + m a J (1 — a)(2 — 3a) 

where we used that p/q = 2/(l + m a ). Note that 2 - 3a - a 2 > and 2 -3a + a 2 - 2a 3 > OwheneverO < a < 1/2, 
so that we must choose e such that only ( |4.29t and A4.3U hold. It is therefore sufficient to let 

1 f a 2 (l-2a)(2-3a-a 2 )l 
£ « = 2 mm { 2-3^' (l- a)(2-3a) j " (4 ' 33) 

It can be easily checked that for any a G (0, 1/2) we have e a > 0. With this choice of e = e a we may take the I 1 
norm of (I4.281 l and obtain that 

9 G Z 2 (J; ^ ( + 1 !T„ iQ)/2 , 1 ) C L 2 (J; ^ + V, )/2ll )- < 4 ' 34 ) 



The Besov embedding theorem C B M / gives that 
so that choosing p = p(a, d) > 2 to satisfy 



Dl+e a Al+e Q -2d/(p+pm Q ) 

jy p(l+m«)/2.1 C ^oo,l » ( 4 - 35 ) 



(4.36) 



p(l + m a ) 
we obtain 

V0eL 2 (X;^ a ). (4.37) 
We note that we may explicitly solve for p in ( 14.361 ) 

V = 7T~r T > 4d > 2 (4.38) 

e Q (l + m a ) 

for any a G (0, 1/2). We recall from ( |4~26l l that V(9 G L 2 (X; L 2 ), and hence V(9 G L 2 (l; 1? n B^i), by d4~37l i. 
Lastly, we use L 2 n 1 C i?^, 1 and the borderline Sobolev embedding theorem 

£W C (4.39) 

to obtain that 

V0 G L 2 (1;L°°) (4.40) 

which concludes the proof of the lemma. Note that by d4.351 l we may even obtain that u G L 2 (2; L°°). □ 

A simple consequence of this improved regularity is the following statement. 
Proposition 4. Let 9 be a Holder continuous weak solution of (11.1 b — ( TT~~4-b such that 

V9eL 2 ([t Q ,t 1 };L°°). (4.41) 

Then 

9 G L°°([t2, ti\; H m ) (4.42) 

for any m > 2, and a.e. i 2 G (*o>ii)- 
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Proof of Proposition^ Since u is divergence free, we have the a priori estimate 



I|||V0|| 2 2 + ||A0|| 2 2 < 



d k Ujd k 6dj( 



(4.43) 



We absorb the (1/2) || A9 1 1| 2 term on the left side of the above estimate and obtain that 

\W)\\%i < II e(t )\\l ie ^ llVe ^ da (4.44) 

which is finite for all to < t G I thanks to the assumption V# G L 2 (X; L°°), as long as 9 (to) £ H 1 . The latter 
is true for a.e. to > since we a priori knew that 6 G L 2 ((0, oo); -ff 1 ), and L 2 functions are finite a.e. (by using 
arguments similar to [;8; Chapter 9], one may even obtain explicit bounds in terms of ||0o||l 2 )- This shows that 
G L co (X; H 1 ) n L 2 (I; H 2 ). Repeating the above argument with X = [to, t{\ replaced by some X' = [t 2 , t\], where 
t 2 > t , we get 

/ "A(u i i 0)A0 



A9\\ 2 2 + \\VA9\\ 2 , < 



<c\\vao\\ L 2\\ao\\ L 2\\vo\\ l ~, (4.45) 



on I' , and therefore obtain 9 G L°°(X';H 2 ) n L 2 (X' ; H 3 ). Hence, for any m > 2, finitely many iterations of the 
above argument proves that 9 G L°°(X"; H m ) n L 2 (X"; H m+1 ) for any X" C X, concluding the proof. □ 

Proof of Theorem\J} If a G (1/2, 1) the theorem follows from the arguments given in Section [3] If a — 1/2, we 
simply consider the solution to lie in C 1 / 2-6 C L°° PI C 1 / 2 , for some e > 0, reducing the proof of the theorem to the 
case a G (0, 1/2). In this case, we apply Lemma[3]to obtain that V9 G L 2 (X; L°°) for any X = [to, ti] C (0, oo). 
Using Propositions] this implies that 9 G L°°([t2, t{\; H m ) for some large enough m (any m > d/2 + 1 is sufficient), 
and a.e. t 2 G (t a , h). The statement of the theorem follows from the Sobolev the embedding H m (R d ) C C l!,5 (R d ), 
for some 6 G (0,1). In particular, if d = 3 one may let m = 3. □ 

Corollary 5. The Holder continuous weak solution 9 of dl.U - dl.4t is C°° smooth for positive time. 

Proof. Let 9 = di9 for some i G {1, . . . ,d}. From Theorem [TJ we obtain that for any < t\ < T we have 
9 G L°°([ti,T];C s (R d )) for some S > 0. From Proposition g] we also have that 9 G L°°([ii, T]; L 2 (R d )) n 
L 2 ([^i, T]; i/ 1 (R d )). Lastly, the equation satisfied by 9, obtained by applying di to dl.lt . is 

d t 9- A9 + i),0'l, ,i),0 ■ i), r,/i;)/> l) (4.46) 

where we use the summation convention over repeated indexes, and are Calderon-Zygmund operators. Given that 
the coefficients d j 9,d t T l] 9 G L°°([t 1 ,T];C 5 (M d )) (cf. Theorem[D are smoother than the a priori smoothness of the 
velocity in ( ll.lt (which belonged to a Holder space of negative index), it is straightforward to repeat the arguments 
used to prove Theorem[TJin order to show that 9 G L°°([t 2 , T]\ C ln (R d )) for some 7 > 0, and a.e. t 2 G (t\,T). 
Since i G {1, . . . ,d} was arbitrary, this shows the solution 9 is C 2n for some 7 > 0. The proof of the corollary is 
concluded by further taking derivatives of the equation, and iterating the above arguments. □ 

5. Higher regularity for the critically dissipative modified SQG equations 

Here we address the applicability of the method presented in Section |4] above, to prove higher regularity for the 
modified critically dissipative SQG equation dl.7t - dl.8t . for the parameter range f3 G (1, 2). Note that when f3 — 2 
the equations dl.7t - dl.8t reduce to the heat equation, and regularity is trivial. In |[T7"1 . Miao and Xue prove the global 
existence of weak solutions 9 G L°°([0, 00); L 2 ) n L 2 ((0, 00); H 13 ^ 2 ), using methods similar to |8][T2]|22], the local 
existence of smooth solutions, and the eventual regularity of the weak solutions (see also 12 [3] @] Q3] and references 
therein, for further results concerning generalizations of the SQG equations). Moreover, in HTl Proposition 5.1], the 
authors prove the following regularity criterion: if a weak solution 9 lies in L°°([to, 00); C a ), with a > (f3 — l)/2, 
then 9 G C°° ((tx, 00) xR 2 )foranyti > to- Such a minimality requirement on a seems to be purely technical, as the 
problem is subcritical in C a for any a > 0. 
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The proof of Theorem[T]of the present paper directly applies to ( |1.7t -( fL"8l , with f3 £ (1,2), and gives the following 
regularity criterion for weak solutions. 

Theorem 6. Let 9q £ L 2 be given, and letT = [to, ii] C (0, oo). Given a weak solution 

9 £ L°°(I;L 2 (M 2 )) nL 2 (l;H p/2 (R 2 )) (5.1) 
of the initial value problem associated to i 1.71 — ( 11.81 1, if 

9 £ L°°(l;C a QSL 2 )) (5.2) 



where 



and (3 £ (1,2), then 



2-/3 (3-1 



<a<l (5.3) 



9£L^([t 2 ,t 1 };C 1 ^(R 2 )) (5.4) 
for some S > 0, andfora.e. t% £ (trj,ii)- Additionally, we have 9 £ C°°((to,ti] X R 2 ). 

Proof. First we note that if a £ ((/? — l)/2, 1), for any /3 £ (1, 2) this result was proven in ifTTI . Therefore, in 
order to complete the proof of the theorem it is left to treat the range (2 — 0) jl < (fi — l)/2, which is equivalent to 
j3 £ (3/2, 2), under the regularity criterion that (2 — /3)/2 < a < (ft — l)/2. In order to avoid redundancy we only 
outline the differences with the proof of TheoremQ] 

Using methods directly corresponding to those described in the proof of Lemma |2] we first prove that if a weak 

b/3/2 
5 P,2 



solution 9 satisfying dBTfl-Ot is such that 9 £ L 2 (l\ B^ 2 2 ), for some p > 2, then 9 £ L 2 (l; Bq/ 2 ), for any 



1 < r < oo, and for any value of q > p such that 



q e( m :\ 4^\- (5.5) 



p \/3 -1- a 13/2 -2a 

Due to our choice of a, the range (15.5b is not empty. Since initially £ L 2 (I; B^)^ we can therefore bootstrap 
this argument finitely many times, and similarly to the proof of Lemma [3] we show that for p large enough we have 
9 £ L 2 (I] B^ 2+2 ^ p ) £ L 2 (I;B^l\). As in Proposition @] this implies, via energy estimates and interpolation 

inequalities, that 9 £ L°°([t 2 , h]; H m ) n L 2 ([t 2l ti]; H m+f} / 2 ) for any m > 1, and a.e. i 2 £ (io,ii). The proof of 
(15.31 follows now from the Sobolev embedding, while the proof of higher regularity consists of taking derivatives of 
the equation and repeating the arguments listed above. □ 

Due to the sub-criticality of the C a norms, for any a > 0, with respect to the natural scaling of the equations 
(11.7b -( [TT8] l, we conjecture that condition ( 15.2b may be replaced with < a < 1, for any j3 £ (1, 2). 

Acknowledgment. The work of S.F. was in part supported by the NSF grant DMS 0803268. The authors would like 
to thank Peter Constantin for fruitful discussions on the topic, and for pointing out the advantage of space-time Besov 
spaces. 

References 

[1] L. Caffarelli and A. Vasseur, Drift diffusion equations with fractional diffusion and the quasi- geostrophic equation. Annals of Mathematics 171 
(2010), No. 3, 1903-1930. 

[2] D. Chae, P. Constantin, and J. Wu, Inviscid models generalizing the 2D Elder and the surface quasi- geostrophic equations. Arch. Ration. 
Mech. Anal., to appear. Available at arXiv:1010.1506vl [math.AP] 

[3] D. Chae, P. Constantin, and J. Wu, Dissipative models generalizing the 2D Navier-Stokes and the surface quasi-geostrophic equations. Avail- 
able at larXiv:1011.0171H [math.AP] 

[4] D. Chae, P. Constantin, D. Cordoba, F. Gancedo, and J. Wu, Generalized surface quasi-geostrophic equations with singular velocities. Avail- 
able at larXiv:1101.3537H [math.AP] 

[5] J.-Y. Chemin, Perfect incompressibe fluids. Oxford Lecture Series in Mathematics and Its Applications, 14. Clarendon Press, Oxford University 
Press, New York, 1998. 

[6] J.-Y. Chemin and N. Lerner, Flot de champs de vecteurs non lipschitziens et equations de Navier-Stokes. J. Differential Equations 121 (1995), 
no. 2,314^328. 



HIGHER REGULARITY FOR PARABOLIC EQUATIONS WITH SINGULAR DRIFT 



11 



[7] Q. Chen, C. Miao, and Z. Zhang, A new Bernstein's inequality and the 2D dissipative quasi-geostrophic equation. Comm. Math. Phys. 271 

(2007) , no. 3, 821-838. 

[8] P. Constantin and C. Foias, Navier— Stokes equations. Chicago Lectures in Mathematics, University of Chicago Press, Chicago, IL, 1988. 
[9] P. Constantin, A. J. Majda, and E. Tabak, Formation of strong fronts in the 2-D quasi-geostrophic thermal active scalar. Nonlinearity 7 (1994), 
no. 6, 1495-1533. 

[10] P. Constantin, G. Iyer, and J. Wu, Global regularity for a modified critical dissipative quasi-geostrophic equation. Indiana Univ. Math. J. 57 

(2008) , no. 6, 2681-2692. 

[11] P. Constantin and J. Wu, Regularity of Holder continuous solutions of the supercritical quasi-geostrophic equation. Ann. Inst. H. Poincare 

Anal. Non Lineaire 25 (2008), no. 6, 1 103-1 1 10. 
[12] S. Friedlander and V. Vicol, Global well-posedness for an advection-diffusion equation arising in magneto-geostrophic dynamics. 

Ann. Inst. H. Poincare Anal. Non Lineaire, to appear. Available at arXiv: 1007. 1211 v2 [math.AP] 
[13] A. Kiselev, F. Nazarov, and A. Volberg, Global well-posedness for the critical 2D dissipative quasi-geostrophic equation. Invent. Math. 167 

(2007), 445-453. 

[14] A. Kiselev and F. Nazarov, A variation on a theme of Caffarelli and Vasseur. Zap. Nauchn. Sem. POMI 370 (2010), 58-72. 
[15] A. Kiselev, Nonlocal maximum principles for active scalars. Available at arXiv: 1009.0542vl [math.AP] 

[16] PG. Lemarie-Rieusset, Recent developments in the Navier— Stokes problem, Reserach Notes in Mathematics, Vol. 431, Chapman & Hall/CRC 
Press, Boca Raton FL, 2002. 

[17] C. Miao and L. Xue, On the regularity of a class of generalized quasi-geostrophic equations. Available at arXiv: 101 1.6214vl [math.AP] 
[18] H.K. Moffatt, Magnetostrophic turbulence and the geodynamo. IUTAM Symposium on Computational Physics and New Perspectives in 

Turbulence, 339-346, IUTAM Bookser., 4, Springer, Dordrecht, 2008. 
[19] F. Planchon, Sur un inegalite de type Poincare. C.R. Acad. Sci. Paris Ser. I Math. 330 (2000), no. 1, 21-23. 
[20] G. Seregin, L. Silvestre, V. Sverak, and A. Zlatos, On divergence-free drifts. Available at arXiv: 1010.6025vl [math.AP] 
[21] L. Silvestre, On the differentiability of the solution to an equation with drift and fractional diffusion. Available at arXiv: 1012.2401 vl [math.AP] 
[22] R. Temam, Navier-Stokes equations. Theory and numerical analysis. Reprint of the 1984 edition. AMS Chelsea Publishing, Providence, RI, 

2001. 

[23] J. Wu, Global solutions of the 2D dissipative QG equation in Besov spaces. SIAM J. Math. Anal. 36 (2004/2005), no. 3, 1014-1030. 

Department of Mathematics, University of Southern California, 3620 S. Vermont Ave., Los Angeles, CA 90089 
E-mail address: susanf ri@usc . edu 

Department of Mathematics, University of Chicago, 5734 University Ave., Chicago, IL 60637 
E-mail address: vicol@math . uchicago . edu 



